The Prediction of a Gapless Topological "Haldane Liquid" Phase in a 
One-Dimensional Cold Polar Molecular Lattice 
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We show that ultracold two-component fermionic dipolar gases in an optical lattice with strong 
two-body on-site loss can be used to realize a tunable effective spin-one model. Fermion number 
conservation provides an unusual constraint that £V (Sf ) 2 is conserved, leading to a novel topological 
liquid phase in one dimension which can be thought of as the gapless analog of the Haldane gapped 
phase of a spin-one Heisenberg chain. The properties of this phase are calculated numerically via 
the infinite time-evolving block decimation method and analytically via a mapping to a one-mode 
Luttinger liquid with hidden spin information. 



I. INTRODUCTION 

We predict theoretically (and support it by numerical 
calculations) the existence of a completely new quantum 
topological gapless liquid phase in certain classes of inter- 
acting one dimensional systems, which, we believe, can 
be realized in ultracold fermionic polar molecular gases 
recently realized in laboratories 1 -. Our predicted novel 
topological phase can in some sense be construed to be 
the gapless analog of the well-known Haldane gapped 
phase that exists in antiferromagnetic spin-one Heisen- 
berg chains^. 

Duan et al£ have recently shown how to realize a 
tunable effective spin-half model with ultracold two- 
component atoms in an optical lattice via superexchange, 
thus allowing, in principle, the experimental study of 
quantum magnetism in the controllable atomic systems, 
but the low temperature scale set by the superexchange 
is currently a barrier to this program. Very recent 
experimental achievement of ultracold gases of dipolar 
molecules 1 - has opened up additional exciting possibili- 
ties for realizing designer Hamiltonians due to the ex- 
tended range of the induced electric dipole-dipole inter- 
action which decays as 1/r 3 . Efforts are currently un- 
derway to load these ultracold molecules into an opti- 
cal lattice. In that case one could realize spin models 
with strong direct off-site interactions, without recourse 
to the inherently weak superexchange mechanism. Our 
proposed new gapless topological phase should be realiz- 
able in such a dipolar molecular fermionic optical lattice. 

One interesting proposal is to use bosonic polar 
molecules to realize an incompressible phase with a hid- 
den order similar to the Haldane gapped phase of a spin- 
one chain 4 . A recent proposal for explicitly constrained 
bosons finds similar results^. This intriguing prospect 
has not yet been experimentally observed, though, as it 
is specific to bosonic molecules while so far only fermionic 
dipolar molecules have successfully been realized near 
degeneracy^. Furthermore, strong two-body loss rates 
for reactive molecules as in Ref£ would severely limit 
the stability of states with multiply occupied sites. 

In this work, we present a mapping of two-component 
dipolar fermions to an effective spin-one model. Rather 



than being an obstacle, the on-site loss rate is actually 
a necessary prerequisite for our scheme. We perform 
numerical calculations for a one-dimensional (ID) lat- 
tice and observe a phase exhibiting hidden order and 
topological degeneracy. By topological degeneracy, we 
mean a dependence of the ground state degeneracy on 
the topology such as occurs in the Haldane gapped phase, 
where the ground state is four-fold degenerate for open 
boundary conditions and unique for periodic boundary 
conditions^. In contrast to the "Haldane insulator" of 
the bosonic proposal^, the phase found here is gapless 
and compressible, and we refer to it as a "Haldane liq- 
uid." 

The layout of the paper is as follows. In Sec. HH we 
write the effective spin-one Hamiltonian. In Sec. IIIII we 
show how the spin-dependent interactions can be experi- 
mentally tuned. Ground state properties of the Hamilto- 
nian are presented in Sec. IIVI with numerical results in 
IIVAI followed by an analytical treatment in llVBl Experi- 
mental detection schemes are briefly discussed in Sec. |V1 
and we summarize in Sec. IVI1 The Appendix contains 
details of the analysis of Sec. IIVBI 



II. EFFECTIVE SPIN-ONE DESCRIPTION 

The two fermionic components (i.e., the two spin com- 
ponents) can be provided by two hyperfine states of 
40 K 87 Rb molecules 8 . The strong on-site loss 6 will give 
rise to effective hardcore repulsion via a continuous quan- 
tum Zeno effect which stabilizes the gas in the presence 
of the lattice^. Thus, after loading the molecules into the 
lattice, the state will quickly decay to a stable configura- 
tion with each site either empty or singly occupied. For 
simplicity, we consider geometries such that the interac- 
tions are isotropic. This is not an essential requirement. 
The low-energy Hamiltonian is then 
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where a =|, I denotes the pseudospin state and Ps is 
the projector onto the subspace with at most one par- 
ticle per site, i.e., J2a c l<j c i^ — 1- Note that in ID 
the hardcore constraint in conjunction with number con- 
servation breaks the Hilbert space up into disconnected 
sectors since a given sequence of t and I fermions cannot 
evolve to a different sequence, although the locations of 
the empty sites may change. 

Since the on-site Hilbert space is spanned by three 
states (t, i, empty), we pursue a mapping to a spin-one 
description, as in Refs. 10 - 11 . To construct spm-one oper- 
ators which commute properly from the physical fermion 
operators we introduce a slave fermion representing the 



empty state, c^o- Then one can verify that 

&i = c tt c »,t~ c !,4. c ^' &t — C l,f C i,0+ c i,O c i,i> $i = (^t) 

(2) 

are effective spin-one operators. 

In terms of the slave fermion, the hardcore constraint 
becomes J2 n c ir) C i,v ~ 1' where rj =t,|, 0, and the 
hopping term in Eq. (fT]) changes as — tc\ a Ci+i, a — > 

— tc\ a c\ +1 Ci+i t(7 Ci t Q. The constraint can be automati- 
cally satisfied by rewriting the Hamiltonian in the form 
of a spin lattice: 
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where the S^-dependence of the hopping ensures conser- states, 



vation of J2i (*- e - 
and annihilation). 



prevents fermion pair creation 



III. ENGINEERING SPIN-DEPENDENT 
INTERACTIONS 



We now show how the interactions can be tuned. In 
the following we keep only nearest neighbor interactions 
and take the lattice to be one-dimensional, but this can 
easily be generalized to further interactions and two di- 
mensions. The Hamiltonian ((3|) becomes especially sim- 
ple for Vff = Vu = —Vfi — V, when it reduces to the 
constrained t — J z model of Ref. 12 for a stripe segment 
in a high temperature superconductor. For 40 K 87 Rb 
molecules in a typical uniform external magnetic field 
the resonant microwave transition frequencies between 
the ground and first excited rotational manifolds are dif- 
ferent for different hyperfine species^. (We will consider 
the Zeeman shift of the noninteracting energy levels to be 
absorbed into the chemical potentials.) Thus one could 
apply an AC field at the average excitation frequency as 
shown in Fig. l(a)| so that the detuning is equal and op- 
posite for the two hyperfine species. If we begin with all 
molecules in the rotational ground state and the field is 
turned on adiabatically, in the limit of large intermolecu- 
lar separation, a molecule will reside in the dressed state 
adiabatically connected to its ground state. The differ- 
ence in the sign of the detuning for the two hyperfine 
species connects their ground states to different dressed 



|er) =e l "*/ 2 cos6»|0, ( 7)-s (T e- 4 " i/2 sin6l|I,a), (4) 

where tan 28 = SI/ 8, SI is the Rabi frequency, 8 is the 
detuning, lo is the field frequency, | J, a) is the bare state 
in the rotational manifold J and the hyperfine state er, 
and Sfj = +1 (—1) for a =| (4). We suppress the mj 
label since the field will dominantly couple the J = 
state to a single J = 1 state whose value of raj is fixed 
by the polarization. 

The resulting effective dipole moments are 
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where E is the polarization unit vector, <i ff = 
dsind cos 9, and d is the transition dipole moment of 
the molecule between the ground and first excited state. 
The time-averaged intermolecular interaction induced by 
a field linearly polarized at an angle 8e to the z-axis is 
then given by 
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More generally, one could control V aa i by carefully tun- 
ing a weak DC electric field and the external magnetic 
field such that the Atoj = transition for one hyperfine 
species is nearly resonant with the Amj = 1 transition 
for the other species, and there are no other transitions 



near resonance, as shown in Fig. 1(b) Then, applying 
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FIG. 1: Sketch of applied AC fields yielding different inter- 
species and intraspecies interactions for two hyperfine states 
in a magnetic field, (a) Scheme yielding Vff = V44 = — Vfj., 
(b) Scheme allowing more general control. 



two microwave AC fields with the same frequency tuned 
near resonance with the chosen transition, one circularly 
polarized and one linearly polarized, each species is af- 
fected by only one of the fields and the interactions can 
be tuned individually. We take the circularly polarized 
field to propagate at an angle 9^ relative to the z-axis 
(which is fixed by the direction of the lattice) and the 
linearly polarized field to have its polarization vector at 
a polar angle 9^ relative to the z-axis and azimuthal an- 
gle <j) relative to the plane formed by the z-axis and the 
propagation vector of the circularly polarized field. As 
before, in the limit of large separation, the molecules will 
reside in the dressed states, but now with il replaced by 
Cla, S by S a = E J= x,a - Ej=o,<7 ~ u, and s a by sgn (S a ). 
The effective dipole moments, d CT = (<r|d|er}, are 

= d c ff 1 (x sin 9^ cos <fi + y sin 9± sin <j> + z cos 9]/) cos cot , 

(7) 

d-j- = d c ff^ (x cos 9^ cos cot + y sin cot — z sin 9^ cos cot) . 



where <L 



icft.a = -sgn (5 a ) d^l 1 - n /+^ 
averaged interaction of molecules on the z-axis is then 
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For intermolecular distances less than or on the 
order of rg = (d 2 /H\5'\) 1 ^ , a full coupled-channel 
Born-Oppenheimer calculation becomes necessary, as in 
Refs. 14 . However, typically the optical lattice wavelength 
A ~ 1/xm, which is greater than rs for detunings as small 
as a few kHz, so we do not need to consider the short- 
range structure of the potential. 



IV. GROUND STATE PROPERTIES 

In the remainder of this work we will consider the 
ground state of the system, taking the experimental pa- 
rameters such that (if = i^i = fi, Vff = V4.4, = V and 
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FIG. 2: Phase diagrams of Hamiltonian © in ID with /i^ = 
[il = fi, Vff = Vu = V and Vfi. = — V cos % and keeping only 
nearest neighbor interaction. V/t strictly equal to zero is not 
included. The range of chemical potential shown corresponds 
to filling factors between and 1. 



Vfi = —Vcosx- (We do not discuss the special case 
cos x = — 1) corresponding to the trivial situation of spin- 
independent interactions.) 



A. Numerical results 

We have calculated the ground state properties of the 
original Hamiltonian ([1]) numerically using the infinite 
time-evolving block decimation (iTEBD) method^. The 
phases shown in Fig. [5] are defined by the structure of the 
entanglement spectrum, as we discuss below, and char- 
acterized by the correlation functions shown in Fig. [3] 

For V > the ground state clearly must lie in the 
sector of the Hilbert space containing states like fCU/fOOl- 
• ■ -, i.e., where the hyperfine state of each molecule is 
different from that of the previous molecule, regardless 
of how many empty sites lie in between. This is long 
range order in the string correlator—, 
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with a = z. (We do not consider the special case of 
V = 0, which is separated from the V > phase by 
the closing of the energy gap between the low-lying z- 
string ordered states and states which break this order.) 
When V cos x > i — nf2^, the entanglement spectrum is 
nondegenerate and every lattice site is occupied, so the 
ground state is an antiferromagnetic insulator (AFM), 
and must be doubly degenerate. In a system with a fixed 
number of molecules less than the number of lattice sites, 
this will be a phase separated region. 

When V cos x < t— M/2 ; there are vacancies and the en- 
tire entanglement spectrum is exactly doubly degenerate, 
which is an indication of a nontrivial topological state**. 
However, in contrast to the "Haldane insulator" of Ref. 4 , 
here the suppression of pair creation/annihilation pro- 
cesses has removed the gap to the longitudinal magnon 
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1. Effective spvnless Luttinger liquid description 
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Consider a lattice with sites 0,1,..., N — 1 and peri- 
odic boundary conditions loaded with n < N fermions 
at sites r = {n < ?'2 < •■• < r n } bearing spins 
(7 = {(Ti, (72, ■■■&n}- Restricting ourselves to the phys- 
ically relevant string ordered subspace, eliminating the 
spin information leaves a spinless attractive fermion sys- 
tem. Any eigenstate of the spinful Hamiltonian, H = 
T + V, in this subspace can be written as a superposi- 
tion *(r,cr) = E 7 =t,; a 7^7 ( r ' °") = E 7 a 7M cr )</'( r )' 

where the two terms correspond to the two possible real- 
izations of string order, i.e., choice of the first fermion's 
spin. This choice uniquely determines all the spins, so 
the state separates into spin and density parts, and 4> ( r ) 
is an eigenstate of the spinless attractive fermion Hamil- 
tonian, H = T + V. (Note that the interaction in the 
string ordered subspace is independent of the two spin 
states.) 



FIG. 3: Correlation functions calculated numerically with 
cosx = 1 for V/t = Q.l,n/t = -1.6 (HL), and V/t = 
1.5, n/t = -0.8 (AFM). 



mode, as can be surmised from the power law decay 
of the spin-spin correlations along z shown in Fig. [3l 
while the transverse modes remain gapped^. Further- 
more, the string correlator along x no longer displays 
true long range order, but decays with a power law. Al- 
gebraic decay of string correlation has previously been 
found in the Hubbard model with infinite onsite repul- 
sion (i.e., Luttinger liquid with Heisenberg exchange) 20 . 
In the present case, though, both long range and quasi- 
long range string order exist simultaneously. Most impor- 
tantly, we shall show below that the doubly degenerate 
entanglement spectrum is the direct result of topological 
degeneracy. We refer to this phase as a "Haldane liquid" 
(HL). 



B. Analytical results 

The mapping from hard-core spin-half fermions to a 
spin-one model provides a suggestive and intuitive lan- 
guage to describe the system, but it is difficult to explain 
the behavior of the HL phase in this language. However, 
given an underlying string- type spin order, one can al- 
ternatively make a mapping from the original hardcore 
spin-half fermions to a spinless Luttinger liquid with a 
prescription for how to access the spin information from 
the density, similar to the approach in Ref . 12 ' 20 . This 
provides analytical solutions for the low-energy physics. 
We heuristically present the mapping below. A more de- 
tailed presentation is given in the Appendix. 



If the lattice is packed, the two string ordered sectors 
are disconnected and all eigenstates are doubly degener- 
ate. For n < N, though, the hopping term, T, connects 
states tpt ( r : f ) an( i ip± (r, a) by taking a particle from 
site N — 1 to site (or vice versa) via the periodic (or 
more generally, phase twisted) boundary conditions. In 
this case the symmetry of the Hamiltonian requires the 
eigenstates to have definite parity, p = ±1, under "time 
reversal" (TR), i.e., T^4> so the eigenstates are unique. 

Consistency of the spinful and spinless Schrodinger 
equations requires TvP (r, a) — E 7 a 7 £ 7 (c) T<fr (r). This 
is clearly satisfied for hops that do not move a particle 
between sites N — 1 and 0. Consistency for hops between 
sites N — 1 and requires e l * = pe 2 *, where <!> (3>) is the 
boundary phase twist of the spinful (effective spinless) 
system. So for TR even (odd) states, we may describe 
the spinful state in terms of an effective spinless state un- 
der boundary conditions with the same (an additional it) 
phase twist. The meaning of this is readily seen in that 
N hops are sufficient to cycle the density around the ring 
back into its original configuration, but the resulting spin 
part of the wavefunction is time-reversed. Thus, TR odd 
states acquire a ir phase after one circuit. One interesting 
ramification is that, if an effective flux is introduced, the 
ground state energy as a function of flux has half the pe- 
riod of a spinless system, as illustrated in Fig. @] by exact 
diagonalization in the string ordered Hilbert space. 

Given the above prescription to obtain the string or- 
dered spinful wavefunction from a spinless one, one ob- 
tains (see Appendix for details) algebraic decay of the x 
string, z spin-spin, and density-density (deviation from 
mean) correlation functions in the HL phase from Lut- 
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FIG. 4: Energy vs. flux for a 5-site ring with 4 string ordered 
fcrmions. Solid (dashed) lines are TR even (odd) states. 



tinger liquid theory: 
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where K > 1 is the Luttinger constant. We have checked 
that the numerical calculations give consistent results 
for K when comparing the various correlation functions. 
Also, the exponential decay of (S§S^) — J2j( c o,-y c r,-f) 1S 
readily understood since one cannot propagate one par- 
ticle past another without destroying the string order, so 
the correlation function must exponentially decay with 
a characteristic length set by the average interparticle 
spacing. 



2. Entanglement spectrum degeneracy 

From the above arguments it is clear that under open 
boundary conditions, the ground state must always be 
doubly degenerate, as also shown in RefJ^. However, for 
periodic boundary conditions it may be unique. To link 
the double degeneracy of the entanglement spectrum to 
this topological degeneracy, consider the reduced density 
matrix for a block of sites, A, obtained by tracing over 
the remaining subsystem, B, 



PA ( r Aj ca; r A> °a) = 



rB,fB 



* (rA,r B ;o'A,crB)1'* (r A ,rB;o- A ,o- B ) ■ (15) 



Since a given number of particles and realization of string 
order in B specifies the number and realization of string 
order in A, the above decomposes as 

/u(rA,cr A ;r A)( 7 A )= © p^ A ' l] (r A , r A ) . (16) 

"A, 7 

In the nondegenerate case, the states have definite TR 
parity, so p^ A '^ = p^ and the entanglement spec- 
trum is exactly doubly degenerate. 



V. EXPERIMENTAL DETECTION 



To probe the HL phase experimentally, one should ini- 
tialize the system in the z string ordered sector (e.g., 
by starting with a 3D lattice and adiabatically turn- 
ing on a strong cylindrical trap). Local correlators 
are directly observable. For instance, (SqS^) can be 
probed by measuring the physical pseudospin-dependent 
single-fermion correlators, J2cr( c l o c i,°)-> v ^ a pseudospin- 
dependent timc-of-flight measurement. Likewise. (S[]S~) 
can be probed by measuring the density-density corre- 
lator, ((/^ (0) - pi (0)) Of (r) - Pi (r))}, via pseudospin- 
dependent time-of-flight noise correlations or repeated in 
situ imaging with single-site resolution. Although finite 
temperature introduces an exponentially decaying enve- 
lope to the gapless correlation functions, at sufficiently 
low temperature the algebraic decay is dominant up to 
some thermal length. 

The entanglement spectrum and string correlators are 
not as directly accessible. However, one can observe other 
consequences of the topological degeneracy. For example, 
in principle, when realizing the simplest case of x — 
with the arrangement of Fig. 1(a) one could use a ro- 
tating ring-shaped optical latticed The half periodicity 



of the ground state energy as a function of rotation fre- 
quency (see Fig. 2]) compared to the corresponding sys- 
tem of spinless fermions could then serve as a signature 
of the phase. 



VI. SUMMARY 

We have shown that an optical lattice loaded with two- 
component dipolar fermions can realize an effective spin- 
one model. Our proposal allows for - in fact, requires 
- experimentally realistic conditions of strong two-body 
loss. In one dimension, we have provided numerical and 
analytical calculations to show that this model contains 
a novel gapless topological phase. 

This work is supported by AFOSR-MURI, DARPA- 
QUEST, and ARO-DARPA-OLE. 



Appendix A 

1. Hamiltonian and Hilbert space decomposition 

Consider a set of spin-1/2 Ising hard core fermions on 
a ID lattice with anti-ferromagnetic nearest neighbor in- 
teractions. (This model has previously appeared in a dif- 
ferent context^.) Let us suppose the ID lattice is a finite 
size ring with some number of sites N. The Hamiltonian 
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may be written as 

H= £ \ C t+l.A l ~ n r+l)Cr,a + tl.C - ^cl a C r . a ] 
r,a— ±1 

- U Cz,aa°'z,l3pCr+'y,a C r,p c r,P c r+-y,a (Al) 

Q,/3,7=±l 

where n(r) — c\ a c r . a . Furthermore let us suppose 
that the chemical potential /i is tuned such that the sys- 
tem has an even number (n) fermions in the ground state. 

Each particle in this system can be labeled by a com- 
posite variable y = (r, a) where < r < N is the po- 
sition of the particle and a = ±1 is the spin state. For 
our purposes we will need to define an order between 2 
values yi — (ri,ai) and y 2 = (^oa) based on the or- 
der of n and r 2 i.e. y\ < y 2 = r± < r%. The Hilbert 
space H of the system above is spanned by the ket states 
U = [|0 < ?/i < y 2 < ■ ■ ■ < y n < (iV-1)}]. The hard core 
constraint together with the Ising constraint on the spins 
leads to an important conserved quantity in the Hamil- 
tonian given by the sequence S = ({crj}j—i f ... t n)(modR) 
i.e the spin sequence is conserved modulo cyclic permu- 
tations R. Thus the Hilbert space may be decomposed 
into a direct sum 

% = ®s^ E (A2) 

where £ represents various equivalence classes of se- 
quences of spins of the n atoms. Each Hilbert space 
7is for a given spin configuration has a simple struc- 
ture Ht. = [|0 < ri < r 2 < ■ ■ ■ < r„_i < r n < N)} such 
that the spin degree of freedom is now implicit. The 
spin-degree of freedom can be restored given the value of 
S. 



2. Low energy wave-functions 

After removing the spin-degree of freedom, each 
Hilbert space and corresponding Hamiltonian have the 
form of spin-less hard core particles moving on a ID 
lattice. The interactions between these particles is at- 
tractive if the spin configuration is in the equivalence 
class £ = (+, —,+,—,...) i.e the state is antiferromag- 
netic. The ground state is therefore purely antiferromag- 
netic. There are 2 antiferromagnetic spin sequences in 
the equivalence class namely, E + = (+, —,+,—,...) and 
= (-,+,-,+,...). 

Thus a wave-function in TLy: can be decomposed into 
2 parts * = + such that * a e Hs a - To un- 
derstand the condition under which "J is an eigenstate, 
we apply the Hamiltonian in Eq. (|A1[) to vj/. For conve- 
nience we consider all position variables lumped together 
as r = (ri, . . . ,r„). All terms in the Hamiltonians un- 
der consideration are diagonal (and therefore like a po- 
tential energy term) other than the hopping term. The 
composite vector r is taken to be in a bounded region 
such that ri < fj+i- The hopping terms T for most 



of the points allow only hopping inside this bounded 
region. Because of periodic boundary conditions, the 
hopping T also connects the boundary points with r 
such that r = (0 < r\ < • • • < r„ = N — 1) and 
r = (0 = ri < • • • < r„ < N — 1). We refer to these hy- 
perplanes as S- and 5+ respectively in parameter space. 
The Schrodinger equation takes the abstract form 

T*(y) = (E - C/(y))*(y) (A3) 

where y = Q < yi < y 2 < ■ ■ ■ < y n < (N - 1), T are 
the allowed nearest neighbor hopping operators in the 
multi-dimensional space and E is the energy eigenvalue. 

Since the spin ordering of each state is alternating, 
the potential energy is independent of the spin part i.e. 
U(y) = U(r). Although the kinetic energy can mix the 
two spin sequences, it is also symmetric, so for n < N the 
eigenstates are nondegenerate and have definite parity 
under time-reversal. Therefore, for n < N, given a state 
^ G Us with time-reversal parity p — ±1, we can define 
a spin-less wave-function ip(r) defined by 

ip(r) = *(0 < (n,+) < (r 2 , -)<■••< (r n , -) < N) 
+ P^(0< (n,-) < (r 2 ,+) < ••• < (r m +) < N) (A4) 

which satisfies the Schrodinger equation 

7>(r) = (E- f/(r))V>(r) (A5) 

where T is the spin-less hopping operator on the lattice. 
The above Schrodinger equation arises from the spin-less 
version of the original Hamiltonian in Eq. (|A1|) . 

H f = y^[4 +1 c r + h.c - A*4c r ] - U 4 +7 4c r c r+7 , 

r 7=±1 

(A6) 

which is the Hamiltonian for spin-less fermions on a lat- 
tice. Here we also choose \x such that the ground state 
has n fermions. Note that the hard core constraint is 
automatically enforced. Similar to our original model 
the Hilbert space can be described by wave-functions 
"0(0 < n < r 2 < ■ ■ ■ < r n < N - 1). 

Conversely, the spin-ful wave-function can be obtained 
given the spinless one ip by 

*(y) = E CaXc»<Mr) (A7) 

a=±l 

where £ + = and \a is the spin part of the wave- 
function and is given by 

Xa(v)= II Sfo + «(-!)')• (A8) 

3=1 n 

The separable form of the wave-function is also obtained 
for the large-[/ limit of the Hubbard model by Refs^^ 
and used by Ref^ in a similar program. Inserting 
Eq. (|A7[) into the right-hand side of Eq. (IA3|) . we obtain 

T*(y) = £ CaXa^f^r). (A9) 

a=±l 
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The validity of Eq. (IA9|) is easy to see in the interior 
of the space r, since the hopping term does not transfer 
any particles across the boundary at r = 0. For points y 
containing particles at the boundaries r = 0orr = iV — 1, 
the equation becomes 

J2 C-aXa(aKV(r) = Yl CaXa(<0e i9 V(r), ( Al °) 

a=±l a=±l 

where 8 (9) is the phase twist at the boundary of the 
spin-ful (spin-less) lattice and ip'(r) denotes the density 
configuration after the hopping operator is applied. So 
we see that for VP time-reversal even, we may describe 
the system in terms of a spin-less wave-function by using 
Eqs. (|A7I) and (| A5|) . where the phase twist on the bound- 
ary of the spin-less Hilbert space is the same as that of 
the spin-ful Hilbert space. But for ^ time-reversal odd, 
the boundary conditions used in Eq. (|A5I) to obtain the 
spin-less wave-function must have an additional n phase 
twist. 

The above arguments are specifically for n < N. For 
the filled lattice n = N, the hopping term does not act 
on the wave-function, and the eigenstates are two-fold 
degenerate. 

3. String order 

The Hilbert space Ha of the low-energy anti- 
ferromagnetic states can be compactly characterized by 
defining a spin-1 operator S r at each site r of the lattice 



such that S r . z = <jj for filled sites rj and SV,z =0 for all 
other sites. With this definition, the Hilbert space Ha is 
characterized by the expectation value of the string order 
of S z operators defined as 

0* tr (r < /) = S r , z J] j* s *<*S r ,, z . (All) 

r<p<r' 

The operator 0* tI (r < r') is non-zero only if both r 
and r' are occupied. The string term in the middle 
Ilr<p<r' e" rSp ' iS computes the parity of the number of 
occupied sites in-between r and r'. Thus 0| tr (r < r') = 
n(r)n(r') where n(r) is the occupation operator at site r. 
The asymptotic expectation value 

(O z stI (r < r')) = (n(r)n(r')) ~ (n(r)) 2 (A12) 

for large r' — r. 

4. Entanglement spectrum degeneracy 

The entanglement spectrum of a block A of length I 
starting from r = 0tor = Z — lof the ring. The rest of the 
ring we refer to as B. In the spin-1 operator notation, the 
Hilbert space is specified by the states |{<SV,2}r=o,...,A r -i) 
where SV, Z are the spin-1 operators at position r. The 
ground state wave-function can be correspondingly writ- 
ten as ^o({S r ,z}r=o,...,N-i)- The reduced density matrix 
is defined as 



PA({Sr,z}r=0, {S' rtZ }r=0,...,l) = ®({Sr,z}r=0,...,N -l)** ({S' r , z }r=0,...,l{S r ,z}r=l+l,...,N -l) ■ (A13) 

{S r ,z}r=l + 1,...,N-1 

I 



The spin variable S r ^ z = for empty sites. Therefore the 
spin-configurations of non-zero weight must satisfy 

£ Si , = S t + E S l* =N-n. (A14) 

r 0<r<l Kr<N-l 
I 



Thus J2o< r <i Sr.z = Eo<r<l &r,z = n A an d the reduced 
density matrix is diagonal in the number of fermion ua 
in the block A. Therefore the reduced density matrix can 
be written in the fermion notation as 



PA(0<yi<---<Vn A <l;0<y' 1 <---<y' nA <l)= *(0<vi<--- <Vn<N) 

l<Vn A +l<—<Vn<N 

x **(0 < y[ < ■ ■ • < y' nA < I < y nA+1 < ■ ■ ■ < y n < N). (A15) 
I 

From the above it follows that pa{® < Vi < ■ • • < y nA ^ respect to time-reversal and thus 

< y'x < ■ ■ ■ < y' nA < I) oc 5 ao .a' - Furthermore, for the 
case with fermion vacancies, \& has definite parity with Pa{® < J/i < ■ • • < Vn A < I', < y\ < • ■ ■ < y' <l) = 

Pa(0 <n< - < r nA < I; < r'j < • • • < r' nA < 1)5^. 

(A16) 
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It follows that the reduced density matrix in the sector 
with riA particles is a direct sum of 2 identical spin- 
sectors. This proves that the entanglement spectrum 
must be doubly degenerate. 



5. Luttinger correlation functions 



b. S z correlation functions 



The model discussed can also be thought of as a spin-1 
Hamiltonian. In this language, the natural correlation 
functions are the S r . z S r ^ z correlation functions. The op- 
erator S r 2 is defined as 



Given that the wave-function of our system can be re- 
lated to the wave-function of spin-less fermions on a lat- 
tice with attractive interactions (Eq. (|A7jl ). we expect 
to be able to relate correlation functions to the spin- 
less fermion correlation functions which are known from 
bosonization of the Luttinger model. The correlation 
functions of operators that preserve the antiferromag- 
netic spin structure (i.e do not introduce spin defects) 
can be mapped to correlation functions of the spin-less 
Luttinger model and therefore turn out to be gapless. 



a. Density-density correlation function 

The density correlation function is the simplest to con- 
nect to the spin-less model. The density operator is de- 
fined by p(r) = S(r — rj).The density-density corre- 
lator can be written in terms of wave-functions as 

(P(r)p(r')) = £ l*o(y)| E S ( r ~ r »W r ' - r ") 

y a,b 

= ^ |^o(r)| X) 5(r - r a )5(r' - r 6 ). (A17) 

r a,b 

The above expression is precisely the density-density cor- 
relation function of spinless fermions on a lattice and is 
given by 



2 A' 



(A18) 

where po = (p(0)) and K > 1 is the Luttinger constant 
for spin-less fermions with attractive interactions^. Here 
A is taken to be the lattice scale. 



(A19) 



The correlator is calculated by evaluating 



y a,b 

= E W)\Y, e l7r(a ~ b) Hr - r a W - r b ) 

r a<b 

= (p{r)e"^-<^' p{p) p(r')} (A20) 



where the last correlator is for the spin-less fermions. 
Here for the second to last statement we have used the 
fact that the spin configuration is antiferromagnetic and 
therefore the spins are opposite if there are an even num- 
ber of particles between a and b. Following ReL2i, we 
use 



p(x) = po--d x <j)(x)+^— 



in the string to get 



(A21) 



(S r ,zSr',z) = {p(r)e-« r 'y^- kF ( r '- r ))p(r')) (A22) 



where for the small lattice constant limit kpa <C 1 we 
have ignored the back-scattering term. Here we also have 
to be careful to take only the real part of this expression. 
Substituting Eq. (IA21[) into the above to express every- 
thing in terms of <fi leads to the expression 



(SrM = p§e- <fc *'*<e iW(r) -* (0)) ) + ( Po + — ^ 



-i(0(r)-<«O))\ 



1 



Po 



-i3kpr 



J3(<f>(r)-<f>(0)) 



-ikpr 



2 e - lkF r ^{r)- m)) + ( po + 



2?rA J 2tt\ n ' V 2yrA / 27rA 

{dr^{r)e lWr) ~ m) dr<p{Q)) - P ° e {(d r (p(r) + d r <t>{0)) e^t'-WW)) 

7T 

ikpr / -y 



i 



2?rA / 2ttA 



Po 



\ e 13 ^ l c i3(<t>(r)-<j,(0))\ 



2ttA J 2ttA 



-ikpr 



+ 



2p e 



■a r (e ,tfW ^ (0)) ). (A23) 



The slowest decay occurs from the first two terms and 
we get 



P 2 o + ^+ ' 



2ttA (2ttA) z 



X\ K/2 
cos (kpr) ( — J 

(A24) 



c. S x string correlator 

In a previous section, we found that the ground state 
was string ordered in the S z string correlations. A nat- 
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ural question to ask is whether this system is also string 
ordered along S x similar to the Haldane model. To as- 
sess this we study the correlation function of the string 
operator along x 



Ol t (r<r') = S r 



n 

r<p<r' 



p ' x S r 



(A25) 



The spin-1 matrix S x is not diagonal in spin space. In- 
stead it converts filled sites to empty sites in the fermion 
representation and empty sites are converted into the su- 
perposition |+) + |— }. Since S x changes the number of 
particles, the action of 0^ tr on is non-zero only if ex- 
actly one of the sites r and r' is filled. The factor e 47rSx 
apart from an overall — sign, flips all spins and leaves 
empty sites unchanged. Therefore the action of Og tr on 
an antiferromagnetic ground state adds a particle 
at r or r' and removes a particle from the other site, 
while flipping the spins in the middle to retain the anti- 
ferromagnetic order. Therefore in the spin-less fermion 
language the string correlator can be written as 

<O s a tr (r < r')> = {-l) {r - r,) (J{r)e™^< v < r > p(p) c (/)). 

(A26) 

Here one must keep the real part only. As before ignoring 
back scattering terms in the exponent 

(0* r (r < r')> = (c t (r) e - l (^( r ')-^W"^('-'- r )) c ( r ')). 

(A27) 



From Ref. 24 , c(r) = e cos (tf>(r) — kpr). Using this 
expression, 



,-t(a^(0)-<^(0)+e(r)-fl(0)) 



s=±l 



+ ( 1 )'' c »2fc F r ^ (| e -«(20(r)+ S 0(O)-0(O)+e(r)-e(O))^ 

(A28) 



s=±l 



Evaluating the expectation values, the leading power-law 
behavior comes from the first term with s = +1: 

(_-t\r / \\ 1 / 2K 
(O* tI (0<r)) = ^-(-\ . (A29) 



6. Gapped correlation functions 

From our discussion of Og tr , the product S rtX S r \ x 
clearly is non-zero only when it acts on configurations 
where exactly one of r and r' is empty. Furthermore, ap- 
plying S r . x S r i. x to the antiferromagnetically aligned state 
^ introduces a defect at the filled site (since it empties 
the site) unless all sites between r and r' are empty. The 
probability of this decays exponentially with an exponent 
proportional to po\r — r'\. The fermion correlation func- 
tion in the spinful fermion model is similarly localized. 
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